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I. INTRODUCTION 



\Q , Extending a method developed by Sasaki in the Schwarzschild case and by Shibata, Sasaki, 

Q\ • Tagoshi, and Tanaka in the Kerr case, we calculate the post-Newtonian expansion of the gravitational 

wave luminosities from a test particle in circular orbit around a rotating black hole up to 0(v 8 ) 
beyond the quadrupole formula. The orbit of a test particle is restricted on the equatorial plane. 
We find that spin dependent terms appear in each post-Newtonian order, and that at 0(v 6 ) they have 
a significant effect on the orbital phase evolution of coalescing compact binaries. By comparing the 
post-Newtonian formula of the luminosity with numerical results we find that, for 30M < r < 100A/, 
the spin dependent terms at 0(v 6 ) and 0(v 7 ) improve the accuracy of the post-Newtonian formula 
significantly, but those at 0(v 8 ) do not improve. 

PACS numbers : 04.25.Nx, 04.30.-w, 04.30. Db, 04.70.-s 
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Among the possible sources of gravitational waves, coalescing compact binaries are considered to be the most 
promising candidates for detection by near-future, ground based laser interferometric detectors such as LIGO [Q, 
VIRGO GEO600, TAMA and AIGO. There are two reasons for this: first, we can expect sufficiently large 
amplitude of gravitational waves from these systems. Second, the estimated event rate, for neutron star binaries, is 
several/yr within 200Mpc ||. Furthermore, the observations of coalescing compact binaries are potentially important 
$H . because they bring us new physical and astronomical information. They can be used to test general relativity [Q, 
J^P" to measure cosmological parameters || and neutron star radii. It may even be possible to obtain information about 
the equation of state of neutron stars If a neutron star or a small black hole spirals into a massive black hole 
with mass < 300AfQ, the inspiral wave form will be detected by the above detectors. Such wave forms carry detailed 
information about the spacetime geometry around the black hole, and therefore may be used to test the black hole 
■ no hair theorem 0. 

When a gravitational wave signal is detected, matched filtering will be used to extract the binary's parameters 
(i.e. masses, spins, ets.) ||. In this method, the parameters are determined by cross-correlating the noisy signal 
from the detectors with theoretical templates. If the signal and the templates lose phase with each other by one 
cycle over ~ 10 3 — 10 4 cycles as the waves sweep through the LIGO/VIRGO band, their cross correlation will be 
significantly reduced. This means that we need to construct theoretical templates which are accurate to better than 
one cycle during entire sweep through the LIGO/VIRGO band Q. If we have accurate templates, we can, in principle, 
determine the mass of the systems within 1% error Thus, much effort has been expended to construct accurate 
theoretical templates M . 

The standard method to calculate inspiraling wave forms from coalescing binaries is the post-Newtonian expansion 
of the Einstein equations, in which the orbital velocity v of the binaries are assumed to be small compared to the speed 
of light. Since, for coalescing binaries, the orbital velocity is not so small when the frequency of gravitational waves is 
in LIGO/VIRGO band, it is necessary to carry the post-Newtonian expansion up to extremely high order in v. A post- 
Newtonian wave generation formalism which can handle the high order calculation has been developed by Blanchet, 
Damour and Iyer jl(|[ll]]. Based on this formalism, calculations have been carried out up to post 5 / 2 -Newtonian order, 
or 0(v 5 ) beyond the leading order quadrupole formula[12-20]. An another formalism is also developed up to 0(v ) 
by Will and Wiseman |l(||2(J which is based on Epstein- Wagoner's formalism pl| , p2[ . 

Although the post-Newtonian calculation technique will be developed and applied to the higher order calculation, 
it will become more difficult and complicated. Thus, it would be very helpful if we could have another reliable method 
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to calculate the higher order post-Newtonian corrections. Recently the post-Newtonian expansion based on black 
hole perturbation formalism is developed. In this analysis, one considers gravitational waves from a particle of mass 
H orbiting a black hole of mass M when fj, <C M. Although this method is restricted to the case when /i <C M, one 
can calculate very high order post-Newtonian corrections to gravitational waves using a relatively simple algorithm in 
contrast with the standard post-Newtonian analysis. This direction of research was first done analytically by Poisson 
p3| who worked to 0(v 3 ) and numerically by Cutler et.al. |2J] to 0(v 5 ). Subsequently, a highly accurate numerical 
calculation was carried out by Tagoshi and Nakamura |25| to 0(v s ) in which they found the appearance of log v terms 
in the energy flux at 0(v 6 ) and at 0(v 8 ). They also clarified that the accuracy of the energy flux to at least 0(v 6 ) 
is needed to construct template wave forms for coalescing binaries. Tagoshi and Sasaki [Ea, using the formulation 



built up by Sasaki 27 1, performed analytic calculations which confirmed the result of Tagoshi and Nakamura. These 
calculations were extended to a rotating black hole case by Shibata, Sasaki, Tagoshi and Tanaka(SSTT) ||8j to 0(v 5 ). 
They calculated gravitational waves from a particle in circular orbit with small inclination from the equatorial plane 
to see the effect of spin at high post-Newtonian orders. They found that the effect of spin on the orbital phase is 
important at 0(v 5 ) order when one of the stars is a rapidly rotating neutron star with its pulse period less than 2 ms 
or a rapidly rotating black hole with q = Jbh/M 2 > 0.2. This analysis was extended to the case of slightly eccentric 
orbits by Tagoshi B9 |, The absorption of gravitational waves into the black hole horizon, appearing at 0(v & ), was 
also calculated by Poisson and Sasaki in the case when a test particle is in a circular orbit around a Schwarzschild 
black hole @. 

In this paper, we extend these analyses in the rotating black hole case to 0(v s ) order. Once again, the calculation 
is based on the formalism developed by Sasaki |?7| to treat a Schwarzschild black hole. Based on the post-Newtonian 
expansion of the luminosity in the test particle limit when the central body is a Schwarzschild black hole ( ]25| J2r| ) , 
Cutler and Flanagan |3l| estimated that we will have to calculate post-Newtonian expansion of gravitational wave 
luminosity at least up to 0(v 6 ) in order to obtain the theoretical templates which cause less systematic errors than 
statistical errors for the LIGO detector. Further, in a previous paper |2j|, we suggested that the effect of spin at 
0(v 6 ) to the orbital phase of coalescing binaries wouldn't be negligible if spin of the black hole was large (i.e. \q\ ~ 1). 
Also, the perturbation study can provide accurate templates for binaries with M 3> /i(for example, binaries of lOOAfo 
black hole-1.4M Q neutron star). Since LIGO and VIRGO will be able to detect gravitational wave signals from 
binaries with masses less than ~ 3OOA/0, it is important to construct templates for such binaries. The frequency of 
gravitational waves from such a massive binary, however, comes into the frequency band for LIGO and VIRGO at 
r/M ~ 16(100M o /M) 2 / 3 , i.e., highly relativistic region. We do not know whether the convergence property of the 
post-Newtonian approximation is good or not in such a highly relativistic motion. Hence, it is an urgent problem 
to clarify at what point the convergence property of the post-Newtonian expansion is good. For these purposes, we 
study the effect of spin beyond 0(v 6 ) order in this paper. 

The paper is organized as follows. In section 2, we present the basic formalism to perform the post-Newtonian 
expansion in our perturbative approach. First we perform the post-Newtonian expansion of the Teukolsky radial 
function using the Sasaki-Nakamura equation. We also show the post-Newtonian expansion of the angular equation, 
which is given in Appendix F. In section 3, we first describe the post-Newtonian expansion of the source terms. We 
consider circular orbits in the equatorial plane around a Kerr black hole. Then the gravitational wave luminosities to 
0(v 8 ) beyond the quadrupole formula are derived. In section 4, we compare post-Newtonian formulas with numerical 
data which gives the exact value of gravitational wave luminosity and investigate the convergence property of the 
post-Newtonian expansion. Section 5 is devoted to summary and discussion. 

Throughout this paper we use the units of c = G = 1. 

II. GENERAL FORMULATION 

A. The Teukolsky equation 

We consider the case when a test particle of mass /i travels in a circular orbit around a Kerr black hole of mass 
M \i. We follow the notation used by SSTT pM), but for definiteness, we recapitulate necessary formulas and 
definitions. 

To calculate gravitational radiation from a particle orbiting a Kerr black hole, we start with the Teukolsky equation 
p2|J||. We focus on the radiation going out to infinity described by the fourth Newman-Penrose quantity, tpi Q, 
which may be expressed as 

/ m e imtp 
duoe- 1 ^ V — - 2 SZ(Q)RemUr), (2.1) 

t.rn ^ 27T 
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where -iS^ is the spheroidal harmonic function of spin weight s = —2, which is normalized as 

' ir |_ 2 5g£| ! Wd0 = l. (2.2) 



The radial function R[ m ui{f) obeys the Teukolsky equation with spin weight s = —2, 

A2 Tr (a~^) ~ ^ Ti muJ (r), (2.3) 

where Ti muJ (r) is the source term whose explicit form will be shown later, and A = r 2 — 2 Air + a 2 . The potential 
Vir) is given by 

K 2 +Ai(r - M)K o . . 

V» = v — — + 8iwr + A, (2.4) 

where K = (r 2 + a 2 )ui — ma and A is the eigenvalue of _2<S^. 

The solution of the Teukolsky equation at infinity (r->oo) is expressed as 

p C„i _^ r 3 e iur f rlirt , Tf, m u{r')R™ muJ {r') _ ~ „3„iwr* / r\ 

ttimuir) -> . in / A2 , n = %»f e , (2.5) 

2 ^ B tmu Jr + A 2 (r') 



where r + = M + V Af 2 — a 2 denotes the radius of the event horizon and -R^L^ is the homogeneous solution which 
satisfies the ingoing-wave boundary condition at horizon, 

ff in f LW,A 2 e _lfcr * for r* -> -oo /„ fi s 

\r- 3 B°^e Mr * +r- 1 S^e- iwr * forr*->+oo, V j 

where fc = cj — ma/2Mr + and r* is the tortoise coordinate defined by 

dr* _ r 2 +a 2 

rfr _ A ' [Z -° 
For definiteness, we fix the integration constant such that r* is given explicitly by 

dr* , 
— — ar 

ar 

rH i-ln i In , (2.8) 

r + - r_ 2M r + - r_ 2M ' v ; 



where r± = M ± V M 2 - a 2 . 



B. Post-Newtonian expansion of the homogeneous solution 

In the previous papers p7| , p8| , the post-Newtonian expansion of the homogeneous solution was performed to 0(e 2 ) 
in the Schwarzschild case and O(e) in the Kerr case, where e = 2Mu>. In this section, we extend those methods, 
performing the expansion of homogeneous solutions up to 0(e 2 ). 

In order to calculate gravitational waves emitted to infinity from a particle in a circular orbit, we need to know 
the explicit form of the source term T^ mw (r), which has support only at r = ro where ro is orbital radius in the 
Boyer-Lindquist coordinate, the ingoing- wave Teukolsky function R^ nuJ (r) at r = ro, and its incident amplitude B 1 ^^ 
at infinity. We consider the expansion of these quantities in terms of a small parameter, v 2 = M/r$. In addition, 
we need to expand those quantity in terms of e = 2Mw since uj ~ 0(Q) where f2 is the orbital angular velocity of 
the particle and Alio ~ 0(v 3 ). In the case of a Kerr black hole, other combination of parameters aui appears in the 
Teukolsky equation. We define q = a/M and we have auo — qe/2 ~ 0(e 3 )- 

First we perform the expansion of the spheroidal harmonics -2Sf£ and their eigenvalues A in terms of aio. Since 
auo = 0(v 3 ), we have to calculate -zSf^ and A up to 0{{au)) 2 ). The eigenvalue A has already been evaluated up to 
0((au}) 2 ) in a previous paper | p8[ . We calculate the expansion of -2<5 , ° I " at 0((auj) 2 ) in the Appendix F. As a result, 
the spheroidal harmonics -2S%£ are given by 



3 



. 2 SZ = -2Pim + <n>S { & + (aufS™ + 0((ac) 3 ), 



(2.9) 



where -iPi m are the spherical harmonics of spin weight s — — 2 fl35fl and 



(2.10) 



Here ci are non-zero only for £' = I ± 1, explicitly 



^im is s iven b y 



^+1 



(£+l) 2 

2 



(£ + 3)(£ - !)(£ + m + !)(£ - m + 1) 



£ 2 



(2£+ l)(2£ + 3) 
+ 2)(^-2)(£ + m)(£-m)]V2 



1/2 



(2€ + 1)(2£ - 1) 



I'm i 



(2.11) 



where the non-zero components of di m are given by 



for any £, by 



d e+1 



d t+2 



rf - - 
a im ~ 



324V7 
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^((O' + lcL 1 ) 2 ), 

(3-m) 1/2 (3 + m) 1 / 2 , 
(3 - m) 1 / 2 (3 + m) 1 /2(4 - m) 1 /2(4 + to) 1 / 2 , 



(2.12) 



for £ = 2, and by 



u £m 



d t+2 

U £m 



120V21 



(4-m) 1 /2(4 + m) 1 / 2 , 



L=(4 - m) 1 /2(4 + m) 1 /2(5 - m) 1 / 2 (5 + to) 1 / 2 , 

180\/n 



it— i 
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(3-m) 1/2 (3 + TO) 1/2 , 



for £ = 3. We don't need for £ = 4 in this paper. 
The eigenvalue A is given by 

A = A + aoj\i + a 2 uj 2 X 2 + 0((aui) 3 ) 
where A = (£ - l)(£ + 2), Ai = -2m(£ 2 + I + A) /{I 2 + £) and 



A 2 



2(£ + i)(c^y + 2e(4-^ + 



2 2 (£ + 4)0?- 3)(£ 2 +£- 3m 2 ) 

3 ~ 3 £(£+ l)(2£ + 3){2£- 1) 



(2.13) 



(2.14) 



Next we calculate the homogeneous solution R™ nuJ - Here we only consider the case when id > 0. We must treat 
the case to < separately. The Teukolsky equation is transformed into the Sasaki-Nakamura equation [^6[, which is 
given by 



r d 2 , . d 
d^- F ^~ U ^ 



X 



£muj 



o. 



(2.15) 



The explicit forms of F(r) and U (r) are given in the Appendix A. The relation between Rimu and Xi muJ is 



4 



Rf.muj — 



(2.16) 



where Ximu — A^ mw A/(r 2 + a 2 ) 1 / 2 , and the functions a, and r\ are shown in Appendix A. Conversely, we can 
express Xi rnu in terms of Remu as 



(r 2 + a 2 ) 1 ' 2 r 2 J_J_ 



(2.17) 



where J_ — (d/dr) — i(K/ A). Then the asymptotic behavior of the ingoing-wave solution Xf^^ which corresponds 
to Eq.(§l| is 



(mi) 1 — i/cr* 



for 



(2.18) 



The coefficient A'^ iuJ , A°^^ and C^ mw are respectively related to B 1 ^^, B^ u and Dg muJ , defined in Eq.(|2.q), by 



B 



pout 



Dtr, 



\ 4™ 

4^ 2 .out 
CO 

1 



(2.19) 



where cq is given in Eq.(A.3) of Appendix A and 



dime = y / 2Mr+[(8 - 2UMu - 16MV)r+ 

+ (12iam - 16M + WamMu + 2AiM 2 uo)r + - Aa 2 m 2 - UiamM + 8M 2 }. 



Now we introduce the variable 



and 



z = z + e 



z+ — Z- 

u>r* + e In e , 



■ In (2 — z + ) ln(z — zJ) 



z+ — Z- 



(2.20) 



where z± — u>r±. To solve A)" by expanding it in terms of e, we set 



X, 



\/ z 2 + a 2 u 2 & m (z) exp (-i<f)(z)) . 



where 



(,)=/*■(§-*) 



* e 1 , z - ^4 
z —z mq in- 



2 z + — Z- z — Z- 



(2.21) 



(2.22) 



which generalizes the phase function to(r* — r) of the Schwarzschild case. This prescription makes it easy to implement 
the ingoing-wav e bou ndary co nditio n on Xf^^. 

Inserting Eq.( [2.2l| ) into Eq.( 2.15 ) and expanding it in powers of e = 2Moj, we obtain 



i (0) fcm] = e£ (1) [& m ] + cQ (1) [&m] + e 2 Q (2) [6m] + e 3 Q (3) [& m ] + e 4 Q (4) fem] + 0(e 5 ), 
where -Z^ 1 -*, Q^ 1 ' and are differential operators given by 



(2.23) 
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dz 2 

z dz 2 
iq\i d 



2 d 

z dz 



£(1 + 1) 



1 


2i\ 


d _ 


('- 


i 


z 2 




dz 


U 3 


z 2 



Aimq 



2z 2 dz 1(1 + l)z 3 ' 



(2.24) 
(2.25) 
(2.26) 



andQ< 2 \ Q( 3 ) andQW are given in Appendix C. Note that the real part of vanishes when we insert the expression 
for Ai. There are A3 or A4 in the formulas for and Q^ 4 \ However, it is straghtforward to show that both A3 and 
A4 do not influence the results in this paper. 



By expanding ^ m in terms of e as 



n=0 



we obtain from Eq.( 2.23 ) the iterative equations 





= 0, 
























+ Q (2) [?0 




















+ Q (2) K&] 





w, 



(3) 



Ira ' 



W, 



(4) 



(2.27) 



(2.28) 
(2.29) 
(2.30) 
(2.31) 
(2.32) 



The general solution to Eq.( 2.281) is immediately obtained as 



AO) _ AO). , AO) 
Qm - a e 3t + Pi 1 



(2.33) 



where jg and rig are the usual spherical Bessel functions. As we discuss later, the boundary condition for n < 2 is 
that ^} is regular at z = 0. Hence (3^ — and we set <xp = 1 for convenience. 

To calculate Q"- for n > 1, we rewrite Eqs.(2.29) — (2.32) in the indefinite integral form by using the spherical 
Bcsscl functions as 



An) 



dzz 2 H W^- U 



dzz 2 n e W^ (71=1,2). 



(2.34) 



The calculation for n = 1 was done in a previous paper j 

Ji)_ (i), (l-W + 3) 
Sim-<*t ]t+ 2{£+l)(2£+l) U 



+ z' 



and we have 

I 2 -A 



2£(2£ + 1 



■Jt-i 



(ndo ~ 3ifi )jo + X] Q + 7^77 ) z2 ( n dk - J**n k )j 



+ n((Ci(2z) — 7 — ln2z) — ^Si(2z) + ijghiz 



imq ( I 2 + 4 

~2~ v^(27TT) 



(l+l) 2 + 4 
(£+l) 2 (2£+l) 



3e+i, 



(2.35) 



where Ci(x) = — dtcost/t and Si(x) — J Q dtsmt/t are cosine and sine integral functions, 7 is the Euler constant 

and a» is an integration constant which represents the arbitrariness of the normalization of Xf^^ 
for simplicity. 

Next we consider Q^j. From Eqs.( 2.34 ) and ( 2. 35] ), and by using formulas in the paper |27), we obtain Q^j as 



We set al 13 = 



f(2) 



.,(2) 



(2), 



„(2) , 



(2.36) 



6 



where /j 2 ' and gf^ are the real and imaginary part of Q 2 ] in the Schwarzschild case respectively, k^l(q) exists only 
in Kerr case and a^2 and are arbitrary constants. The explicit forms of fj; 2 ' and g^ are given in a previous 
paper ^7|. The term k^(q) is given for £ = 2 by 

( 2) _ 191 i . m 2 q 2 j () mqji 

68 i q 2 73m 2 q 2 7mqj3 i 2 * 22 

mqio 70 7 2 H 73 H to a 7s 

63 J S92- 7 1764 J 180 72 y ^ 324 y Ji 

Hi q 2 ji 71m 2 q 2 j 4 13 i 

+ 420 Wgj4 "W" 8820 + ~ mqm 

mqji I3mqj 3 \ f -i . 13 z . \ 

In x+ — mqji - -^rmqj 3 b(Z) 



90 J V 5 " J 90 

5 ' 90 



j 13 z 

— mqni + — — mqn 3 I C(z), (2.37) 



and for I = 3 



(2) 3527 i 2 to ji mqj 2 5i 2 5i 2 2 

& 3m = o;n m ? Jl ^75 77T7 1 32 + TO (J J 2 



840 JJ 315 36 504 " J 2268 

379 z q 2 jz 7 to 2 q 2 j3 3mqj4 i 2 i 22 

W W9J3 ~W" ^2^ + ^6^-140 9 J4 + Tl20 m 9 j4 

97 i g 2 j 5 17to 2 9 2 j 5 103 i 25 i 

to q 75 to q nn H mqn-2 

5040 360 5040 48 8 

13mqj 2 5mqj 4 \ f-13i . 5i . \ 

-VX~ + -56~ J hl Z + I W m 9 J2 " 56 m 9 J4 J S(Z) 
13? hi \ 

— mqn 2 + — mqm J C(z), (2.38) 

where C(z) = Ci2z — 7 — ln2z and S(z) = Si2z. Note that to obtain above two formulas, we have added terms 
proportional to ji to simplify the formulas of below. As noted previously, the source term T^ muJ has support 

only at r — r and ujr — 0(v). Hence we only need X^ nuJ at z = 0(v) < 1 to evaluate the source integral, apart 
from the value of the incident amplitude A™ mul . Hence the post- Newtonian expansion of corresponds to the 

expansion not only in terms of e = 0(v 3 ), but also z by assuming e <C z -C 1. In order to evaluate the gravitational 
wave luminosity to 0(v s ) beyond the leading order, we must calculate the series expansion of Q™ in powers of z for 
n = to I = 6, for n = 1 to £ = 5, for n = 2 to I = 4, for n = 3 to £ = 3 and for n = 4 to £ = 2 (See Appendix C of 
SSTT). 

When we evaluate iL,, we examine the asymptotic behavior of at infinity. Since the accuracy of -A^L we 

need is 0(e 2 ), we don't have to calculate and £^ m closed analytic form. We need only the series expansion 

formulas for Q 3 ^ and Q 4 ^ around 2; = 0, which is easily obtained by Eq.( 2.34 ). Inserting £^ into Eq.( 2.21 ) and 
expanding it by z and e assuming f«z«l,we obtain 

( 3 ) -q 2 i 3 -i 7 rag i 2 2 mq 3 to 3 q 3 

K ' = to q H m q H 

S2m 30z 30z * 30 180 60 y 36 90 

to q In z i 2 2! 

■ to q In z 



30 30 

G 



/ 319 100637 i q 2 17 to 2 q 2 83 i , 

+z ! TOq 1 1 mq 

1 6300 441000 180 1134 5880 



61 2 33 In z 106 z 



ra q + — — — to q In z — — to q In z' 



13230 * 15 1575 * 30 

+0(z 2 ) + a 2 3 J lj2 +f3 2 3 J l n 2 , (2.39) 

C£2 = + O^' 1 ) + "2 4 ii2 + /& 2) (2.40) 
.(3) — i to 2 q 2 i 3 i 3 3 

**»=i260 m 9+ lW- 1260 m<Z ~ 3780 W 9 + ° {z) 
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(2.41) 



The boundary condition of that correctly represent the boundary condition of X^^ (Eq.(2.18)) is that z^"^ 



must be no more singular than z^ +1 at 

and n in this paper. As for o^, they still remain arbitrary and we set a£| 

tin) 



(n) (n) 

0. Since we need Q m only up to n = 4, we set = for all of £ 



for all of I, m and n = 3,4. 



Inserting Q ' into Eq.(2.22) and expan ding it in terms of e = 2Mu, we obtain Xf^^ which are shown in Appendix 
C. Using the transformation of Eq.( 2.16| ), we obtain B^ mu} which are also shown in Appendix D. 

Next, we consider A™ mu at 0(e 2 ). Using the relation ~ —je-i ~ {—^) l+a ni n -i at z ~ oo, etc., we obtain the 
asymptotic behavior of Q*] and £^ at z ~ oo as 



:(2) 



(p£1 + til In z - (In z) 2 ) + - p« ^ *)n< 
+*pim^ ln z + - ln z)n t ln z 



where 



(i) 



7T 

2' 



(1) _ 



(£-!)(£ + 3) 



-In2- 



2imq 
£ 2 (£+l) 2 ' 



1 

I 



fc=i 



(2.42) 
(2.43) 

(2.44) 
(2.45) 

(2.46) 



for any ^ and 



(2) _ 



P-2 



o (2) 

<<2m 



457 7 
210 ~ 

-7 ln 2 

-457 7T 
420 



(2 ) = 527 



7^ 5tt 2 
2 + 24 

" 18 

77T 

f T ^^6 
5^ 2 



— 7 to q 
18 ' y 



457 ln 2 
210 



to g ln 2 - 
5 to q 



(ln2) 2 



3G 



mir q 



2 2 i 

m q + 



(2) 

?3m 



-7 ln 2 

-26 7T 



7' 

1 7 TO o + 

2 24 72 ' q 

i , n (ln2) 2 
to g ln 2 



72 
52 ln2 
21 



324 



7T ln2 



72 

7 7r 67 mq 



2 ' 
i i , 

TO 7T (7 H (7 

144 y 360 



17 i 



2 „2 



7T In 2 



■ to g 



21 2 1440 144 3 360 12960 2 
Then noting that exp(— i<f>) ~ exp(— i(z* — z)) at z ~ oo, the asymptotic form of -2Q™ W is expressed as 

= ^ 2 +« 2 ^ 2 expH0) {/£> + e£« + e 2 ^ + -} 
e--* (^f >e«) [l + e (p« + iq^) + ^ + ^ 



(2.47) 
(2.48) 

(2.49) 
(2.50) 



+e^*(z^ 1) e- lz ) 



(2) 



(2.51) 



where h} and /iV^ are the spherical Hankel functions of the first and second kinds, respectively, which are given by 



,(2) 



h) — jt + 



From these equations, noting wr 

1 



.4 



tmuj 



„ , 1 -ieln 

2 



(-l) w — , = J<-in<-(-l) 
elne, we obtain 

W-l.-,,'",^,-,,'-' — l ;,/2» 



£+1 ' 



(2.52) 



(2.53) 



The corresponding incident amplitude B™ mljJ for the Teukolsky function are obtained from Eq. ( 2.1£ ) . 
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III. GRAVITATIONAL WAVE LUMINOSITY TO 0(V 8 ) 
A. The geodesic equations 

In this section,we solve the geodesic equation for circular motion in the equatorial plane. The geodesic equations 
in the Kerr geometry are given by 



(XT 



C - cos 2 e\a 2 (l - E 2 ) + = 6(0), 

sin > - 



= _ (aE - a sin 2 9 + ( E(r 2 + a 2 ) - al z ) = T, 

dr \ sin 0/ A V / 

Z^=±VR, (3.1) 

(IT 

where E, l z and C are the energy, the z-component of the angular momentum and the Carter constant of a test 
particle, respectively. S = r 2 + a 2 cos 2 and 

R = [E(r 2 + a 2 ) - al z } 2 - A[(Ea - l z ) 2 + r 2 + C]. (3.2) 

Since we consider a motion of a particle in the equatorial plane 8 = n/2, we can set C = 0. We define the orbital 
radius as r = r . Then E and l z are determined by R(r ) — and dR/dr \ r =r = as 

1 - 2v 2 + qv 3 

E = 



(l-3w 2 + 2gw 3 ) 1 /2 ' 
r v(l - 2qv 3 + q 2 v 4 ) 
(1 - 3v 2 + 2qv 3 y/ 2 



where v = (M/ro) 1 ^ 2 - After these preparations, we can easily obtain tp(t) as 

ip(t) = nt , 



M 1 / 2 

n = [1 - qv 3 + q 2 V 6 + 0(v 9 )] . (3.3) 



B. Integration of the source term 

Using results of the previous section, we can now derive the source term of the Teukolsky equation and integrate it 
to give the amplitude of the Teukolsky function at infinity. 

The energy momentum tensor of a test particle of mass n is given by 

T " v = ^dmri d -i s{r - ^ - ^ - «'»■ < 3 - 4) 

The source term of the Teukolsky equation is given by 

/QdLO 
dfldtp-Zp- 1 ^ + B '*)e- imv+wt ^j^, (3.5) 

where 



2tt 
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B' = 



- 5 p 8 pL_l[p- 4 Lo(p- 2 p- 1 rnn)] 

P 8 pA 2 L_ 1 [p- 4 p 2 J + (p- 2 p- 2 A- 1 T l7m )], 



2\/2 



- - P 8 pA 2 J+ [p- 4 J+ (p^pT^)] 

p 8 p^J + \p- i p 2 H- 1 L- l {p- 2 p- 2 T inn )], 



2V2 



(3.6) 



with 



p = (r — ia cos 9) 
vn 

L s = dg H — : — - — auj sin 9 + s cot 6, 
sm9 

J+ = d r + iK/A, 



(3.7) 



and p denotes the complex conjugate of p. 

In the present case, the tetrad components of the energy momentum tensor, T nn , T^n and Tfnm, take the form, 



sin 

Cm 



T- - — 



Kill 

sin 9 



l -6(r-r )6(9~Tr/2)6(v-v(t)), 



%-r )^-7r/2)%-^)), 



(3.8) 



where 



4£ 3 i 



Cm n 



MP 



2 r 



^-'■nr in 



P£_ 

2£i 



i +a 2 )-al z \ , 
[£(r 2 + a 2 ) - al z ] 
I, 



i sin 9 1 a_B 



sin 2 9 



i sin ( aE 



sin 2 9 



(3.9) 



and t = dt/dr. 



Substituting Eq.(3.6) into Eq.(3.E) and integrating by parts, we obtain 

4 



Ti 



£muj 



V2W- 



dt / d9e lut - mv{t) 



1 



--Llip-iLl^S^C^p-'p-'Sir - r )S(9 - tt/2) 

A 2 =2 

+ -7^(4-9 + ia(p-p) sm9S)j + {C r nnP~ 2 r 2 A- 1 6(r - r )S(9 - tt/2)} 
V2p 

+ ^Llip^S^p-^^C-nnAp-'p-'Sir - r )S(9 - tt/2) 
- ]p 3 A 2 SJ + {p- 4 J+(pp- 2 C filr nS(r - r )S(9 - tt/2))} 



(3.10) 



where 



Lt =d 



sin 6* 



auj sin# + scot 8, 



(3.H) 
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and S denotes _2<S^(#) f or si mplicity. 
We further rewrite Eg. ( 3. 10 ) as 



T tmu = / dte^-^^A 2 [(^0 + A mn o + A„- ir - n0 )S(r - r ) 
+ {(4 in i + Ajnihi)8(r - r Q )} + {A r f irfl2 S(r - r )} 



=7T/2 



where A nn Q etc. are given in Appendix B. Inserting Eq.( 3.12 ) into Eq.(2.5), we obtain Z^ mu> as 

2it5(tjj — mil) 



Ztmw — 



2iujB l e 



^(mu{^« n + A rn n Q -+- A m m 

0} 



dR 



(irnuj 



in 

dr 

5(uj — mfl)Zi r< 



A r 



2 pin 



d 2 R 



dr 2 



trau) a 



r—ro,d—7r/2 



(3.12) 



(3.13) 



Using characters of -2S^(9) at 9 = n/2, it is straightforward to show that T^_ m) _ w = (— l) £ T£ imjW where T^ TO , W is 
the complex conjugate of T^ m ,u- Since the homogeneous Teukolsky equation is invariant under the complex conjugate 
followed by m — > — m and u> — ► — w, we have Zi t - m ,-u> = {—l) Zt,m,u>- 



C. Results 



I n thi s section, we calculate the gravitational wave luminosity up to 0(v s ) beyond the quadrupole formula. From 
Eq.(2A), ^4 a t r — > °o takes a form, 



6 I 



qCLUJo 



(3.14) 



£=2 ro=- 

where <x>o = 7nf2. At infinity, ■04 is related to the two independent modes of gravitational waves h+ and h x as 

04 = ^(h+ - ih x ). (3.15) 

From Eqs.( 3.13 ), ( 3.14 ) and ( |3.15 ), the gravitational wave luminosity is given by 

I 2 

I Zg mu!o 



dE 
~dt 



— = V 

^ 47rWo ^ 



dE\ 



In order to express the post-Newtonian corrections to the luminosity, we define r\g, m as 

1 ( dE\ 



dE 
~dt 



N 



where (dE/dt)^ is the Newtonian quadrupole luminosity: 

32/i 2 M 3 32 /^\2 



dE 
~dt 



N 



5rX 



5 \MJ 



We only show r)^ m for m > mode since f]^ m = rji 

107 v 2 



/4784 



?72,2 = 1- 1 ^— + (4 7r-6g) ^+(^: + 2 (? 2 U 4 



21 

-428 7T 4216 q 
21 + 189 



V1323 

/99210071 1712 7 16tt 2 
V 1091475 lOfT + 3 



(3.16) 



(3.17) 
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8830 g 2 3424 In 2 1712 In v\ , 

- 28 nq + v 6 

q 567 105 105 ' 

(191 T7 , 163928(7 t . , ; , > 



+ 



27956920577 183184 7 1712tt 2 207167rg 
81265275 + 2205 63 + 189 

456028 q 2 4 366368 In 2 183184 In 

9261 +q + 2205 + 2205 
v 2 qv*( 17 <?\ 4 

(* 793g\ 5 / 2215 Trg 859g 2 , (i 
\18 9072/ W V 254016 6 1512 1 

-17 7T 11861 g 7Tg 2 7q 3 \ 7 



252 190512 8 12 

/ 15707221 107 7 tt 2 1045 irq 
^26195400 ~ 945 + 27 4536 
118943 q 2 q 4 107 In 2 107 In v 



+ 



190512 16 945 945 
Putting together the above results, we obtain (dE/dt)e = ^2 m (dE/dt)i, n for I = 2 as 

dE\ (dE\ f 1 1277 v 2 , 73g\ 3 , ^37915 , 33 



dt 7 Wi/jvl 252 V 12/ V 10584 16 



1 TTTTT- + M 7T - — H « + TTT^TT + 



+ 



2561tt 201575g\ 5 /2116278473 1712 7 16 



126 9072 7 \ 23284800 105 



„6 



169 irq 73217 g 2 3424 In 2 1712 Inv 

6 + 4536 105 105 ; 

/ 76187 tt 376387 q 65 Trg 2 151 g 3 \ 7 

+ V 5292 + 27216 + 8 12 / W 

2455920939443 548803 7 5129 tt 2 70877 irq 
+ t^ttttt^ + „„„ 77777" + 



v 8 



7151344200 6615 189 648 

64835431 g 2 17 g 4 219671 In 2 548803 In v 
1333584 + 16 + 1323 + 6615 



For I = 3, we obtain 



_ 1215w 2 1215u 4 /3645tt 1215g 
773,3 ~ 896 112 + V 448 112~ ' V 



5 



/ 243729 3645 g 2 \ 6 /-3645tt 131949 g\ 7 
+ V 9856 + 896 J V + \ 56 + 1792 J V 

( 25037019729 47385 7 3645 tt 2 32805 ?rg 
+ V 125565440 1568 + 224 448 

346275 g 2 47385 In 2 47385 In 3 47385 In v 
+ 14336 1568 1568 1568 ' 
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_ 5v 4 40gw 5 / 193 80g 2 
773,2 ~ "63 189~ + V 567 + 567 ' ' 



6 



20tt 352g\ 7 / 86111 1607rg 40g 2 
"63" + 1701 J V + V 280665 189 + "27" 



v 8 . 



V 2 V 4 ( TT 17g 

773,1 ~ 8064 ~ 1512 + \ 4032 ~ 9072 ' V 



5 
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/_437_ + iV\ 
V 266112 24192 / 



+ - 



1137077 
50854003200 
13 In 2 13 1nu 



~* , 3601 9\ v 7 
756 435456 / 

7T 2 145 irq 41183 g 2 

42336 + 6048 ~ 36288 + 3483648 



137 



42336 42336 



Then we obtain 



(di) 3 



+ 



+ 



~dt) 



N 



1367 v 2 32567 v 4 ( 16403 tt 896 q 



/152122 



+ 



1008 
341 q 2 



3024 V 2016 



81 



-13991 tt 



V 6237 ' 81 
/ 5712521850527 79963 7 

V 28605376800 2646 
11168371 g 2 79963 In 2 



435456 



2646 



216 

6151 tt 2 
+ 378 " 
47385 In 3 
1568 



4019665 q\ 



54432 
192005 Trg 

2592 
79963 Inv 
2646 



For I = 4, we have 



151808 v 6 ( 10240 tt 12800 q 
6237 + V 567 567 



+ 



1280 v 4 
567 

/ 560069632 5120 g 



+ 



774,3 



^4,2 



774,1 



V 6243237 ' 567 
729 v 6 729 qv 7 



1792 
437 v e 

+ 



+ 



24640 
20 tt 80 q 



4480 
5v 4 

3969 " 43659 V 3969 '- } >- m 
( 7199152 200 g 2 \ 8 



28431 3645 q 2 
+ 



14336 



V 218513295 

v 6 qv 



282240 112896 



27783 / 

7 / 10] 

+ 



+ 



5<7 2 \ 



4656960 903168 



Then we obtain 



dE 
dt 



(dE\ f 8965 v 4 
\dt) N { 3969 3492720 
/ 51619996697 66084895 q 2 



84479081 v 6 ( 23900 tt 59621 q 



V 1323 



V 582702120 7112448 



2592 



For £ = 5 we have 



7?5,5 



V5A 



7?5,3 



775,2 



9765625 v b 

2433024 
4096 v s 



2568359375 v 8 
47443968 



13365 
2187 v 6 
450560 

4 v s 



775,1 = 



40095 

v 6 



150903 v 8 
2928640 



179 v 8 



127733760 2490808320 



Then we have 
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dE\ _(dE\ f 1002569 w 6 3145396841 v i 



dt J, \dt N \ 249480 58378320 



(3.35) 



For f = 6 we have 



26244 v* 

%6 = ^575- ' (3 - 36) 
131072 v s 

V6A = ^555975- ' (3 ' 37) 

% ' 2 = 573^585 ' (3 - 38) 



and 776,5, r/6,3 , ??6,i become 0(f 9 ). Then we have 

'dE\ ( dE\ 210843872 v 8 



dt J 6 \dt J N 28667925 



(3.39) 



Finally, gathering all the above results, the total luminosity up to 0(v ) is expressed as 

dE\ ( dE\ f 1247v 2 / 73q\ , / 44711 33<? 2 \ , 



dt / \ dt J N { 336 V 12 / V 9072 16 

-8191 tt 3749 c?\ 5 /6643739519 1712 7 



672 336 J V 69854400 105 

16 7T 2 169 7to 3419 q 2 3424 In 2 1712 In v\ fi 

H ^ H v b 

3 6 168 105 105 1 

-16285 7T 83819 q 65irq 2 151 q 3 \ , 

1 w < 



+ 



504 1296 8 12 

323105549467 2325977 1369 tt 2 3389 nq 124091 q 2 



3178375200 4410 126 96 9072 



17g 4 39931 In 2 47385 In 3 2325971nu\ R l 

H — H 1 w 8 > . (3.40) 

16 294 1568 4410 J J y ' 

In Appendix G, we present formulas for ?y^ m and dE/dt in terms of v' = (Mfi) 1 / 3 for the sake of convenience to 
calculate the phase function for a inspiraling wave form l a|. 

Setting q = 0, above reproduces the previous results |£5| (2^] in a Schwarzschild case. Up to 0(v 5 ), the results 
agree with those obtained by SSTT [^8| in the case when the test particle moves a circular orbit in the equatorial 
plane. For I — 5 and 6, there are no contributions due to the black hole spin and the results are identical to the 
Schwarzs child case. 

In Eq.(|3.40|), the numerical value of terms at order 0(v 6 ) is given by w 6 (115.7 - 88.48g +20.35g 2 - 16.30 Inv). We 



find that the spin dependent terms are not so small compared to the other two terms if |g| is of order unity. Thus, 
we see that spin dependent terms at 0(v e ) will give a significant effect to template wave forms of coalescing binaries 
when spin of a black hole is large. 

Finally we note that the angular momentum flux can be easily calculated from 



dJ \ _ 1 IdE 

~dt / ~ n\~dt 



(3.41) 



IV. COMPARISON WITH NUMERICAL RESULTS 



As discussed in section I, it is important to investigate the detailed convergence property of the post-Newtonian 
approximation. Therefore we compare the formula for dE/dt, derived above, with numerical results and investigate 
the accuracy of the post-Newtonian expansion of dE/dt. 

In this section, we consider the total mass of the binary systems including black holes ~ 2 — 300M© because 
gravitational waves from such binaries can be detected by LIGO and VIRGO. In particular, we pay attention to the 
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accuracy of post-Newtonian formula for dE/dt when r < 100M (or v > 0.1), because gravitational waves from these 
binary systems will be detected when the orbital separation becomes less than r ~ 100M. Here, we ignore the effect 
of absorption of gravitational waves by the black hole. We will briefly discuss its effect in the next section. 

A numerical study of dE/dt from a particle in a circular orbit in the equatorial plane around a Kerr black hole has 



been performed by Shibata |39|. Since nothing was assumed about the velocity of a test particle, those results are 
correct relativistically in the limit /i <C M. In that work, dE/dt was calculated with accuracy < 1CP 4 . However we 
found that this accuracy is not sufficient to compare it with the post-Newtonian formula for dE/dt including terms 
up to 0(v a ). Thus, in this paper, we calculate dE/dt again requiring the accuracy to be ~ 10 -5 . In the numerical 
calculations, we have taken into account the contribution from the 1 = 2 through i = 6 modes in dE/dt which is 
consistent with the post-Newtonian formula. 

In figs.l(a— e), we show the error in the post-Newtonian formulas as a function of the Boyer-Lindquist coordinate 
radius when q = —0.9, —0.5, 0, 0.5 and 0.9. In these figures, we show the error for 6 < r/M < 100. Since the radius 
of the inner stable circular orbit for q = 0.9 is ri so ~ 2.32M and a stable circular orbit is possible for r > r; so , we also 
show the errors in the case when q = 0.9 for 2.5 < r/M < 12 in fig. 2. The error in the post-Newtonian formula is 
defined as 



Error : 



! dE\ I (dE\ 

dt /pn/ V dt /nr. 



(4.1) 



where (dE/dt)pw and (dE/dt)-^^ denote the post-Newtonian formula and the numerical results respectively. As for 
(dE/dt) PN , we have used 2-PN, 2.5-PN, 3-PN, 3.5-PN and 4-PN formulas. Here, we define n-PN formula as the 
expression for dE/dt which includes post-Newtonian terms up to 0(v 2n ) beyond the quadrupole formula. In each 
figure, open square, filled triangle, open triangle, filled circle, and open circle denote the error of 2-PN, 2.5-PN, 3-PN, 
3.5-PN and 4-PN formulas, respectively. We note that in fig. 2, the errors in the 2.5-, 3- and 4-PN formulas become 
greater than unity for very small radius, because in such a region, dE/dt for those PN formulas becomes negative. 
From these figures, we find the following. 

(1) If we use the 2-PN or 2.5-PN formula, the error is always greater than 10~ 4 when r < 100M irrespective of q. If 
we use the 3-PN formula, however, the error decreases significantly, and it becomes less than 10 for r > 60M , and 
less than 10 -3 for r > 30M irrespective of q. 

(2) If we adopt the 3.5-PN formula, the accuracy becomes better than that of 3-PN formula. The error is always less 
than 10 -4 when r is greater than ~ 30M and less than 10 -5 when r is greater than ~ 60M. This feature does not 
depend on q. However, if we use the 4-PN formula, the accuracy is not improved compared with the 3.5-PN formula. 
In particular, this tendency is remarkable for smaller radius. 

(3) The accuracy of the 3.5-PN or 4-PN formula is not always better than that of the lower-PN one inside r c , where 
r c < 5M for q = 0.5 and 0.9, r c ~ 10M for q = and —0.5, and r c ~ 15M for q = —0.9. Thus, the convergence of 
the post-Newtonian expansion seems rather poor around r c . 

Using the above results, we investigate the accuracy of the post-Newtonian formulas as templates for various binary 
systems. As explained in section 1, to investigate the accuracy of the post-Newtonian formulas as templates, it is 
useful to check if they can predict the number of cycles of the gravitational waves, N, with accuracy less than 1. 
compact binary systems, the cycles are mainly accumulated around ~ 10Hz which is the lowest frequency region in 
the LIGO band, and N is approximately given by 

iV-LgxW^fV^y (4-2) 



V M ) V 4 /v 

where M and (i are the total mass and reduced mass, respectively. This means that the template must have an 
accuracy less than 

^xHT'f^-fV^, (4.3) 

when the frequency of gravitational wave becomes 10Hz. 

First we consider equal mass binary systems, that is M = 4/i. At 10Hz, the orbital separation of a binary of 
total mass M is approximately given by r/M ~ 347(M /M) 2 / 3 . We find that the 2-PN and 2.5-PN formulas are 
insufficient if M > 5M , and the 3-PN formula is needed. The 3-PN formula seems adequate irrespective of q. 

On the other hand, the situation is slightly different in the case when a neutron star of mass ~ 1.4M Q spirals into 
a larger black hole. In such a case, the number of the cycles of the gravitational waves is large compared with the 
equal mass case when the total mass is the same. Thus, it seems that we need at least the 3.5-PN formula for binaries 
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of mass greater than ~ 3OM to obtain the required accuracy. Also, for binaries of mass greater than ~ 70Af Q , we 
need higher post-Newtonian corrections beyond 4-PN order. 

Binary systems of total mass greater than ~ 100M Q can be detected when r is smaller than ~ 15M. However, 
as mentioned in (3) above, the convergence property of the post-Newtonian expansion becomes bad for small orbital 
separations. In particular, for q ~ — 1, the accuracy of the post-Newtonian expansion seems bad at r ~ 15M. Thus, 
it may not be appropriate to use the post-Newtonian approximation for binaries of total mass ~ IOOA/q with large 
mass ratio ji <C M. A more detailed investigations of the convergence of the post-Newtonian expansion will require 
the calculation to be carried beyond 4-PN order. 



V. SUMMARY AND DISCUSSION 



In this paper, we have performed a post-Newtonian expansion of gravitational waves from a particle in a circular 
orbit around a Kerr black hole. The orbit lies in the equatorial plane and the calculations are accurate to 0(v 8 ) 
beyond the quadrupole level. We have performed the post-Newtonian expansion of the Sasaki-Nakamura equation 
and obtained the Green function of the radial Teukolsky equation up to 0(e 2 ) using methods developed previously. 
Then we obtained all the necessary radial functions to the required accuracy. We have also calculated the spin 
weighted spheroidal harmonics up to 0((au>) 2 ). The outgoing wave amplitude of the Teukolsky function and the 
gravitational wave luminosities were derived up to 0(v 8 ) beyond the quadrupole formula. 

It is worth noting that in the formula for 772,2 in Appendix G, there are terms such as (— 8/3)qv' 3 , 2q 2 v' 4 , (—8/3)q 3 v' 7 
and q A v' 8 . In a previous paper j28|, we pointed out that the term 2q 2 v' 4 can be explained in terms of the quadrupole 
formula as the contribution of the quadrupole moment of the Kerr black hole to the orbit of the test particle. A 
similar explanation is possible for (— 8/2>)q 3 v' 7 and q 4 v' 8 . We can derive those terms by using the quadrupole formula 
dE/dt = 32/5/i 2 f 4 f2 6 , where f is the orbital radius of a test particle in dc Donder coordinates. If multipole moments 
of the black hole exist, the orbital radius is changed due to the influence of those multipole moments (or if we fix the 
orbital radius, is changed due to the multipole moments of black hole). We can calculate the leading order effect 
of the multipole moments to the orbital radius by using multipole expansion of the Kerr metric (Eq.(10.6) of Ref. 

In this way, we find that the dominant effect of the multipole moments of a Kerr black hole to dE/dt can be 
expressed as [Ell 



f - 1 (£)> {' - §*■• - + iS *>" + (4* + H ■"} 



(5.1) 



where Mg and Sg are mass and current multipole moments of a Kerr black hole given by Mi + iSt = M(ia) . Now 
we can interpret the term — \2q 3 v 7 as the effect of the current octopole moment of a black hole and the term q 4 v 8 as 
the effect of both the mass quadrupole moment and t = 4 mass multipole moment of a black hole. 

As for £ = 2 and m = 1 mode, there are terms -qv' 3 /12, q 2 v /4 /16, -7q 3 v' 7 /2A and q 4 v' s /16. The terms -qv' 3 /l2 
and q 2 v' 4 /16 can be explained as the correction to the radiative current quadrupole moment Q |4^| . We expect that 
the terms -7q 3 v' r /24 and <?V 8 /16 can also be derived simply in a similar way. 

In section 4, by comparing post-Newtonian formulas for dE/dt with numerical data, we indicated that the con- 
vergence of the post-Newtonian expansion seems bad when orbital radii of binaries become less than ~ lbM. This 
suggests that the post-Newtonian expansion may not be appropriate to construct theoretical templates for large mass 
ratio binaries where the total mass is greater than ~ IOOMq because gravitational waves from such binaries enter 
the LIGO/VIRGO frequency band when r < 15M. Nevertheless, the higher order post-Newtonian terms gradually 
improve the accuracy of the templates. Hence, it is very natural to ask whether the post-Newtonian expansion is 
always appropriate or not, and if appropriate, up to what order we need the post-Newtonian terms to construct 
accurate templates. Fortunately, it is possible to obtain the formulas for dE/dt which include post-Newtonian order 
terms beyond 0(v 8 ) by extending techniques developed in this paper. Extension of the present work up to the higher 
post-Newtonian order, beyond 0(v s ), is very important and that is our future work. 

The analysis, in this paper, has been restricted to the case when a test particle moves in a circular orbit on the 
equatorial plane. However, as shown in a previous paper p8|| , inclination of the orbital plane from the equatorial 
plane will significantly affect the orbital phase evolution. Hence, the present work should be considered as a first step 
toward the complete calculation of the energy and angular momentum luminosities including the orbital inclination. 

Finally, we comment on the effect of absorption of gravitational waves by the black hole event horizon which should 
be taken into account when we consider the orbital evolution of black hole binaries. According to Gal'tsov 0], the 
lowest order contribution of the gravitational wave absorption to dE/dt is given by 



dE 
It 



(f)iH'+^H}(^ 2 > 



(5.2) 
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Thus, the effect of absorption appears from 0(i> 5 ) if q 7^ 0. Although the coefficient is small compared with that of 
dE/dt for the outgoing wave even in the case \q\ ~ 1, we need the expression for dE/dt due to the black hole absorption 
to obtain an accurate template up to 0(v 8 ). Therefore, to obtain the higher order post-Newtonian corrections to the 
black hole absorption is a problem for the future. 
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APPENDIX A: THE FORMULAE OF F AND U 



In this Appendix we show the potential functions F and U of the SN equation (2.15). Details of the derivation are 
given in Ref. 

The function F(r) is given by 

nr) = V --^- 2 , (Al) 

where 

77 = c + ci/r + c 2 /r 2 + c 3 /r 3 + c 4 /r 4 , (A2) 

with 

c = -\2iuM + A(A + 2) - I2aut(au! - to), 
Ci= 8ia[3acj — \(au> — to)], 

c 2 = -24iaM(auj - to) + 12a 2 [l - 2(auj - to) 2 ], 
c 3 = 2Aia 3 (auj - to) - 24Afa 2 , 

c 4 = 12a 4 . (A3) 



The function U(r) is given by 



where 



U(r) 



G 



( r 2 +a 2)2 



G 



2 _|_ AG r 



FG, 



(A4) 
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APPENDIX B: FUNCTIONS IN THE SOURCE TERM 
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APPENDIX F: SPHEROIDAL HARMONICS 



In this appendix, we describe the expansion of the spheroidal harmonics -vS®!^ at order O{(aoj) 2 ). 
The spheroidal harmonics of spin weight s = —2 obey the equation, 



1 d 



sin 9 d9 



f sin* 4) 

l del 



9 , , „ (to — 2 cos( 
a u> snr 9 - - n 



sm 



Aauj cos 9 — 2 + 2mau> + X 



-20 lm — U. 



We expand -2Sf£ and A as 
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where _2-P| m are the spherical harmonics of spin weight s = — 2 and A n are given in section 2.2. Here 
represent the ^-dependence of A n for later convenience. We set the normalization of -iPgm as 
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where < ji,j2, mi,ni2\J, M > is a Clebsch-Gordan coefficient. Then, for £ = 2 and 3, we obtain df m W^i) which are 
given in section II. B. As for df m , we consider the normalization of -iPim (Eq.( |2.2| )). Inserting Eq.(F2) into Eq.( ^.2| ), 
and using the orthogonality of we obtain 
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APPENDIX G: THE EXPRESSION OF THE LUMINOSITY BY MEANS OF THE ORBITAL ANGULAR 

FREQUENCY 

For the sake of convenience to calculate the orbital phase error, we describe the formula of gravitational wave 
luminosity by means of v' = (A/51) 1 / 3 . In this appendix, we define r)^ m as 
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Figure Captions 

Figs.l(a-e) : Error of the post-Newtonian formulas as a function of the Boyer-Lindquist coordinate radius r for 
6 < r/M < 100 in the case q = —0.9, —0.5, 0, 0.5 and 0.9. In each figure, open square, filled triangle, open triangle, 
filled circle, and open circle denote the error of 2-PN, 2.5-PN, 3-PN, 3.5-PN and 4-PN formulas, respectively. 
Fig. 2 : Error of the post-Newtonian formula of q = 0.9 for 2.5 < r/M < 12. Open square, filled triangle, open triangle, 
filled circle, and open circle denote the error of 2-PN, 2.5-PN, 3-PN, 3.5-PN and 4-PN formulas, respectively. 
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